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ABSTRACT 

A continuous fiber composite is modelled by a two-element composite 
cylinder in order to predict the elastoplas tic response of the compos- 
ite under a monotonically increasing tensile loading parallel to fib- 
ers. The fibers and matrix are assumed to be elastic -perfectly plastic 
materials obeying Hill's and Tresca's yield criteria, respectively. The 
present paper investigates the composite behavior when the fibers yield 
prior to the matrix. 


INTRODUCTION 

The elastoplas tic response of fibrous composites has been the subject of a 
number of theoretical studies [ 1 -4] . When a composite is subjected to uniaxial 
tension loading parallel to the fibers, a two-element composite cylinder (Figure 
1-a) has been frequently utilized to model the composite response. The loading 
direction together with the axisymmetric geometry of the representative volume 
element simplify the mathematical difficulties associated with the equilibrium 
equations. By implementing a traction- free boundary condition to the outermost 
surface of the representative volume element, it becomes possible to construct a 
well -posed boundary value problem when the fibers and the matrix are assumed to 
obey Hill's and Tresca's yield criteria, respectively. A closed form analytical 
solution requires further simplifications such as elas tic-perfectly plastic con- 
stituents, perfect interfacial bonding, etc. When the composite constituents are 
assumed to obey the modified yield criterion proposed by Hill[l], the hardening 
effect of the matrix can be taken into account without mathematical difficul- 
ties. However, the present study focuses on a composite with non-hardening con- 
stituents . 
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(a) 



Figure 1. Continuous Fiber Composite 


(a) Representative Volume Element 

(b) Cross section 
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Hill proposed a relatively simple yield criterion which assumes the differ- 
ence between the axial stress and the arithmetic mean value of the radial and 
circumferential stresses to be equal to the yield stress[l]. When this yield 
condition is implemented to a composite with elastic fibers surrounded by an 
elas tic -perfectly plastic matrix, the entire matrix yields simultaneously, 
resulting in a discontinuity in the slope of the effective stress - strain curve. 

Mulhern, Roger, and Spencer [2] proposed a rigorous analytical solution for a 
two-element composite cylinder under cyclic loading. Their study models an elas- 
tic core fiber surrounded by an elastic-perfectly plastic matrix tube obeying 
Tresca's yield criterion. The resulting composite behavior is almost as bilinear 
as Hill s solution. However, the slope of the effective s tress - strain curve is 
continuous because the plastically deformed matrix zone propagates from the 
fiber -matrix interface to the outer surface of the matrix. 

Ebert and Gadd[3] studied a similar problem for an elastic-perfectly plastic 
core fiber surrounded by elastic matrix. Ebert, et al.[4] extended this to an 
elastic-perfectly plastic matrix. However, the application of their numerical 
solution is restricted to a composite in which the Poisson's ratios of the fiber 
and the matrix are equivalent. 

The present paper extends the study of Ebert and other authors [3, 4] to a 
two-element composite cylinder representing a transversely isotropic fiber sur- 
rounded by an isotropic matrix in which the Poisson's ratios of the core fiber 
and the matrix need not be identical . 


MODEL FORMULATION 

Consider a metal matrix composite reinforced by continuous fibers under uni- 
axial tension loading parallel to the fibers. The globally averaged stress state 
of the representative volume element is assumed to be one dimensional. The elas- 
tic-plastic response of the bar can be analytically predicted by solving an 
equivalent boundary value problem of a single core fiber which is perfectly 
bonded to the surrounding matrix tube. The volume element representing the 
equivalent boundary value problem is illustrated in Figure 1. The uniaxial ten- 
sion loading in the fiber direction produces a three dimensional stress state in 
both the core fiber and the surrounding matrix. When the tension loading 
increases monotonically , either the fiber or matrix yields at a certain magni- 
tude of the applied tension loading. Further increment of the tension loading 
induces the yielding everywhere in the composite. The possible yield sequences 
for the composite constituents of the representative volume element can be cate- 
gorized into three cases as shown in Figure 2. The present study provides ana- 
lytical solutions to the first case under monotonically increasing tension load- 
ing. 



Figure 2. Possible Composite Yield Sequences 
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The entire mathematical formulation of the present study is based on the 
following key assumptions: 

1. The core is assumed to be a transversely isotropic fiber surrounded by an 
isotropic matrix. 

2. Both constituents are assumed to be elas tic-perf ectly plastic. 

3. The interfacial bonding between the core fiber and the matrix is assumed to 
be perfect throughout deformation. 

4. The core fiber is assumed to obey Hill's yield criterion. 

5. The matrix tube is assumed to obey Tresca's yield criterion. 

6. The axial strain is assumed to be spatially homogeneous. 

Since the geometry of the representative volume element is axisymmetric and 
the loading direction is parallel to the core fiber, the only non* trivial 
equilibrium equation is 

do j- <Tj- - o 0 

+ = 0 ( 1 ) 

dr r 


The constitutive equations for the transversely isotropic core fiber and the 
isotropic matrix are 
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respectively . 

Since the plastic strain is incompressible, 


p p p 

e r + €0 + e z = 0 (3) 

It can be shown that Hill's yield criterion becomes the following yield condi- 
tion for the transversely isotropic core fiber under a transversely isotropic 
loading : 



This mathematical expression is identical to Tresca's yield criterion. 


The surrounding matrix will yield according to one of the following conditions 
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external boundary conditions and interior compatibility are 


a r = 0 at r 

- b 


(5-a) 

a r — unique 

at r 

= a or r = a, c 

(5-b) 

u r = 0 at r 

= 0 


(5-c) 

u r = unique 

at r 

= a or r = a, c 

(5-d) 


When both fiber and matrix responses are elastic, the stress state and the 
displacement field in the representative volume element are determined by match- 
ing the radial stress and displacement at the fiber-matrix interface. The prob- 
lem solving procedure for this elastic deformation is simple and straightforward 
as discussed below. 

Since the stress state in the fiber is always transversely isotropic, 
f f 

a r = a e = - p (6-a) 

f 

a z = - 2 "LT p + e L £ z (6-b) 

where P is an unknown constant to be determined. The radial displacement in 
the fiber is given by 

u r = C l r (7-a) 

The strain components become 
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From (2-a),(6), and (7-b), 


Ci = 


(1-j/ tt -2i/ lt i/ tl )P/E t + ^LT e zJ r 


Within the matrix, 


m 


Pa z 


b 2 - a 2 


b 2 

1 - — 

,-2 


(7-b) 


(8-c) 


(9-a) 


58 



m 

On = 


?a z 


m 


b^-a^ 


2P*'m a2 

b 2 !^ 


.2 i 


1 + 


+ E m e z 


u r = C2^ + C3/r 
The strain components become 


(9-b) 

(9-c) 
( 10 - a) 


e r = u r , r = c 2 ‘ c 3/ r 

III r\ 

e 6 = u r/ r = c 2 + c 3 / r 

From (5-a) and (5-b) , C 2 and C 3 are determined as functions of P. 

Pa 2 (l+^ m ) (l-2i/ m ) 
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Pa 2 b 2 ( l+i^ m ) 


( 10 -b) 

( 10 -c) 


( 11 ■ a) 


C 3 = - (H-b) 

E m (b 2 -a 2 ) 

Then the radial displacement in the matrix becomes 
Pa 2 ( 1 +i^ m ) 

u r - — [(l- 2 , m )r + b 2 /r] - u m e z ( 12 ) 

E m (b 2 -a 2 ) 

The interfacial stress in the radial direction, -P, is determined from 
the displacement compatibility given by (5-d) . 
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It can be shown that the effective axial Young's modulus of the volume element 
is given by 

E c = E L (a/b ) 2 + E m (b 2 -a 2 )/b 2 + 2 (a/b) 2 (^ m -^ LT ) 2 /Ai (14) 

Under monotonically increasing axial strain, e z , either the core fiber or the 
surrounding matrix yields first. The yield strains, strengths, and the Poisson's 
ratios of the constituents govern the composite yield sequence. When the applied 
axial strain reaches a certain value, e z y, the core fiber yields first if the 
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longitudinal yield strength of the fiber is much smaller than the yield strength 
of the matrix. The initial yield strain, e zY , at which the entire fiber yields 
is given by 



e L a 1 + ( 1/ m" l/ LT) ( 1_2 i/ Lt) - 


(15) 


After the core fiber yields, the surrounding matrix behaves elastically 
until matrix yield at the interface occurs. Since the plastic strain is incom- 
pressible and the stress state is transversely isotropic in the core fiber, 
stresses, total strains and plastic strains in the fiber become 
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where -P is the unknown interfacial stress to be determined. 

f f f 

The first strain invariant, e r +«^+e z> together with the uniqueness of the 
radial displacement at the material interface determine the magnitudes of the 
fiber stress components as functions of the applied axial strain. During this 
strain increment where the matrix still deforms elastically, the stress compo- 
nents and the radial displacement in the matrix given by equations (9-a), (9-b), 

and (9-c) are still valid if the interfacial stress in the radial direction is 
redefined as 


P = 
where 

n - 


-n I'll 

r?i+2r/ 2 

E_ 


(l' 2t/ m) (1- 2 ^LT) 


^LJ 


(b 2 - a 2 ) /a 2 


l + "m L 1 ‘ 2l/ m + ( b / a ) 2 


E T 

r) 2 ” 

2 - 2j/ tt - ^xl/ v LT 

The radial displacement in the fiber becomes 


1 'mn +, ?2~ 

c _ r + 

(1-2^lt)t72 

y l 
■£ 

rji+2r) 2 

c z 

rii+2r) 2 

e l 


(17-a) 
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Further increase of the applied axial strain causes yielding of the matrix 
material. The plastically deformed region then propagates toward the outer sur- 
face of the matrix. During this strain increment, the elastically deformed mat- 
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rix and the plastically deformed matrix exist together as shown in Figure 1-b. 
Outside the interface between the plastic matrix and elastic matrix, the stress 
state and deformation are given by the elasticity solution with an unknown 
internal pressure, P*. acting on the interfacial surface between the plastic 
matrix and elastic matrix. Therefore, within the elastically deformed matrix, 
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where P' c is determined from Tresca's yield criterion given below. 
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where c is an unknown function of the applied axial strain, e z . 

Within the plastically deformed matrix, Tresca's yield criterion given by 
eqs . (21-a) and (21-b) determines the stress state and displacement field. If 
the Poisson's ratio of the matrix is smaller than that of the fiber, the fiber- 
matrix interfacial stress in the radial direction is always tensile. On the 
other hand, when the Poisson's ratio of the matrix is greater than that of the 
fiber, the fiber-matrix interfacial stress may be either compressive or tensile. 
After the core fiber deforms plastically, the apparent Poisson's ratio of the 
fiber approaches 0.5 as the applied axial strain increases. The interfacial 
stress in the radial direction may be changed from compressive to tensile before 
the initiation of matrix yield. If the matrix yield strain is far greater than 
that of the core fiber, the interfacial stress in the radial direction at the 
onset of matrix yield becomes tensile. Then, from eq. (3) and (4-c) , 
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From eq. (24) and the strain- displacement relationships, the radial and circum- 
ferential stresses are expressed in terms of the radial displacement and its 
gradient with respect to r. 
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The equilibrium equation thus becomes 
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The solution to the above differential equation is given by 
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Since the radial stress and displacement must be single -valued at r=c and 
r=a , the unknown constants are determined as 
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The radius of the matrix yield front, c, can also be determined as a function of 
the applied axial strain, e 2 , and material properties. However, it is more 
convenient to express the applied axial strain as a function of the radius of 
the matrix yield front by satisfying the uniqueness of the radial displacement 
at the fiber-matrix interface. Until the yield front reaches the outermost sur- 
face of the matrix tube, the axial strain is given as 
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After the entire matrix yields, it can be shown that eqns . (28-a) and (28-b) 

should be corrected for further axial strain increment as 
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Then the stress components in the matrix become: 
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Further increase of the axial strain, as mentioned in ref. [2], may cause 
another type of plastically deformed matrix region in which the radial and cir- 
cumferential stresses are identical. The present paper, however, does not con- 
sider this case because the infinitesimal strain assumption may not be valid for 
further increase of the applied axial strain. 
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RESULTS AND DISCUSSION 


The effective stress - strain curve for a composite cylinder can be predicted 
by calculating the average value of the axial stress, a z , as a function of the 
applied axial strain, e z , and the mechanical properties of the composite 
constituents. As an example, the effective stress-strain curve of the composite 
studied by Ebert, et al . [4] is demonstrated in Figure 3. The mechanical proper- 
ties of the composite constituents appear in Table 1. In Figure 3, the solid 
lines represent the present analysis. Within the straight line segment(OA) , both 
the core fiber and the surrounding matrix tube are within their elastic limits. 
When the applied axial strain reaches € zY , the entire core fiber yields. The 
next line segment(AB) represents the hardening behavior of the composite with 
plastically deformed core surrounded by an elastic tube. When the applied axial 
strain reaches e z ^ in the same figure, the surrounding matrix starts yielding 
from the fiber-matrix interface. This strain can be calculated from eq. (29-a) 
by setting c=a. Then the plastically deformed matrix region propagates outward 
until the entire matrix tube yields. This smooth transient region is represented 
by the line segment BC . The applied axial strain, e z2 , where this transient 
phenomenon terminates can be calculated from the same equation by setting c=b . 
The composite response to further axial strain increase then follows the remain- 
ing line segment. Within the transient region and for the higher value of the 
applied axial strain, the matrix tube material is assumed to be nonhardening 
even though the material hardens signif icantly(Figure 3 in ref. 4). The exper- 
imental results of Ebert, et al.[4] are also plotted in the same figure. 


Table 1. Constituent Properties the Composite [4] 


Material 

Ultimate 

Strength 

(Ksi) 

0.05% 

Offset Yield 
(Ksi) 

Elastic 

Modulus 

(Msi) 

Poisson' s 
ratio 

SAE 4140 
(Core ) 

93 

54.9 

28.7 

0.29 

Maraging Steel 
(Matrix Tube) 

318 

288 

25.5 

0.29 


In Figure 4, the radial variations of the radial, circumferential and axial 
stresses in the composite cylinder of which the core volume fraction is 0.5 are 
plotted for two distinctive axial strain values, € zl and e z2 . The stresses in 
the core material decrease slightly as the axial strain increases from e i. At 
the onset of the initiation of the matrix yield, the axial stress in the matrix 
tube is constant. As the applied axial strain increases beyond e z ^, the axial 
stress in the matrix has its maximum value at the free surface. 
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Figure 3. 


Composite 



Figure 4. Stress 







CONCLUSIONS 


The present study provides an analytical prediction of the elastoplastic 
response of continuous fiber composites with weaker fibers. The incremental form 
proposed by Ebert, et al.[4] must be replaced by the second order ordinary dif- 
ferential equation given by eq . (26). Furthermore, the present analysis can 
handle the mismatch of the Poisson's ratios as well as transversely isotropic 
fibers. The present analysis will be generalized for the same type of compos- 
ites under cyclic loading for providing a comparison to the study of Mulhern et 
al. [21 
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